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We study the stability of (4+ l)-dimensional charged Gauss-Bonnet black holes and solitons. We observe 
an instability related to the condensation of a scalar field and construct explicit "hairy" black hole and 
soliton solutions of the full system of coupled field equations. We investigate the cases of a massless scalar 
field as well as that of a tachyonic scalar field. The solitons with scalar hair exist for a particular range of 
the charge and the gauge coupling. This range is such that for intermediate values of the gauge coupling a 
"forbidden band" of charges for the hairy solitons exists. We also discuss the behaviour of the black holes 
with scalar hair when changing the horizon radius and/or the gauge coupling and find that various scenarios 
at the approach of a limiting solution appear. One observation is that hairy Gauss-Bonnet black holes never 
tend to a regular soliton solution in the limit of vanishing horizon radius. We also prove that extremal 
&Q' Gauss-Bonnet black holes can not carry massless or tachyonic scalar hair and show that our solutions tend 

to their planar counterparts for large charges. 

(n: 
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c r i 1 Introduction 

^LJ ' The gravity-gauge theory duality pQ has attracted a lot of attention in the past years. The most famous example 
is the AdS/CFT correspondence [5] which states that a gravity theory in a d-dimensional Anti-de Sitter (AdS) 

7—1 ' space-time is equivalent to a Conformal Field Theory (CFT) on the (d — l)-dimensional boundary of AdS. 
ILJ , Recently, this correspondence has been used to describe so-called holographic superconductors with the help of 
black holes in higher dimensional AdS space-time [3J |U [5) |B] . In most cases (3 + l)-dimensional black holes with 
planar horizons (k — 0) were chosen to account for the fact that high temperature superconductivity is mainly 
^3 ' associated to 2-dimensional layers within the material. The basic idea is that at low temperatures a planar 
black hole in asymptotically AdS becomes unstable to the condensation of a charged scalar field. The hairy 
black hole is the gravity dual of the superconductor. The main point here is that this instability occurs due 
to the fact that the scalar field is charged and its effective mass drops below the Breitenlohner-Freedman (BF) 
bound [7] for sufficiently low temperature of the black hole hence spontaneously breaking the U(l) symmetry. 
Surprisingly, however, the scalar condensation can also occur for uncharged scalar fields in the <i-dimensional 
planar Reissner-Nordstrom-AdS (RNAdS) black hole space-time [3]. This is a new type of instability that is not 
connected to a spontaneous symmetry breaking as in the charged case. Rather it is related to the fact that the 
planar RNAdS black hole possesses an extremal limit with vanishing Hawking temperature and near-horizon 
geometry AdS2 x M d_2 (with d > 4) [HI HI IH] • For scalar field masses larger than the d-dimensional BF bound, 
but smaller than the 2-dimensional BF bound the near-horizon geometry becomes unstable to the formation of 
scalar hair, while the asymptotic AdS^ remains stable [4]. The fact that the near- horizon geometry of extremal 
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black holes is a topological product of two manifolds with constant curvature has led to the development of the 
entropy function formalism [TTJ [T21 [T3J . In [Tl] the question of the condensation of an uncharged scalar held on 
uncharged black holes in (4 + 1) dimensions has been addressed. As a toy model for the rotating case, static 
black holes with hyperbolic horizons (k — —1) were discussed. In contrast to the uncharged, static black holes 
with hat (k = 0) or spherical (k — 1) horizon topology hyperbolic black holes possess an extremal limit with 
near-horizon geometry AdS2 x H 3 and it has been shown by numerical construction that the black holes form 
scalar hair close to extremality. 

Higher curvature corrections appear naturally in the low energy effective action of string theory [15] . In 
more than four dimensions the quadratic correction is often chosen to be the Gauss-Bonnet (GB) term, which 
has the property that the equations of motion are still second order in derivatives of the metric functions. As 
such explicit solutions of the equations of motion arc known. The hrst example of static, spherically symmetric 
and asymptotically flat black hole solutions in GB gravity were given for the uncharged case in [TBI HZ] an d 
for the charged case in [TB]. Moreover, the corresponding solutions in asymptotically Anti-de Sitter (AdS) 
[TJ5J [5D1 US [52] as well as de Sitter (dS) space-times [23j have been studied. In most cases, black holes not 
only with spherical (k — 1), but also with flat (k = 0) and hyperbolic (k ~ —1) horizon topology have been 
considered. Moreover, the thermodynamics of these black holes has been studied in detail (2TJ [22j [24] and the 
question of negative entropy for certain GB black holes in dS and AdS has been discussed [201 ES] . 

Uncharged d-dimensional Gauss Bonnet black holes with hyperbolic horizon topology possess also a regular 
extremal limit such that the near horizon geometry contains an AdS2 factor. For scalar fields with masses 
above the d-dimensional BF bound, but below the 2-dimensional BF bound one would expect the near-horizon 
geometry to become unstable to scalar hair formation. This was demonstrated in [26 . In particular, it was 
shown that the radius of the near-horizon AdS2 decreases with increasing Gauss-Bonnet coupling and tends to 
zero in the Chern-Simons limit [26 . 

Interestingly, there seems to be a contradiction between the holographic superconductor approach and the 
Coleman-Mermin- Wagner theorem [27] which forbids spontaneous symmetry breaking in (2 + 1) dimensions 
at finite temperature. Consequently, it has been suggested that higher curvature corrections and in particular 
GB terms should be included on the gravity side and holographic GB superconductors in (3 + 1) dimensions 
have been studied [28]. However, though the critical temperature gets lowered when including GB terms, 
condensation cannot be suppressed - not even when including backreaction J2H1 [301 EH] • 

In [33] static, spherically symmetric black hole and soliton solutions to Einstein-Maxwell theory coupled to a 
charged, massless scalar field in (4+l)-dimensional AdS space-time have been studied. The existence of solitons 
in global AdS was discovered in [33], where a perturbative approach was taken. In |32j it was shown that 
solitons can have arbitrarily large charge for large enough gauge coupling, while for small gauge coupling the 
solutions exhibit a spiraling behaviour towards a critical solution with finite charge and mass. The stability of 
RNAdS solutions was also studied in this paper. It was found that for small gauge coupling RNAdS black holes 
are never unstable to condensation of a massless, charged scalar held, while for intermediate gauge couplings 
RNAdS black holes become unstable for sufficiently large charge. For large gauge coupling RNAdS black holes 
are unstable to formation of massless scalar hair for all values of the charge. Moreover, it was observed that for 
large gauge coupling and small charges the solutions exist all the way down to vanishing horizon. The limiting 
solutions are the soliton solutions mentioned above. On the other hand for large charge the limiting solution is a 
singular solution with vanishing temperature and finite entropy, which is not a regular extremal black hole [34] . 
These results were extended to a tachyonic scalar field as well as to the rotating case [35]. Recently, solutions 
in asymptotically global AdS in 4 dimensions have been studied in [36]. It was pointed out that the solutions 
tend to their planar counterparts for large charges since in that case the solutions can become comparable in 
size to the AdS radius. 

In this paper, we are interested in the condensation of a charged tachyonic or massless scalar fields on 
charged static black holes and solitons in (4+l)-dimensional Anti-dc Sitter space-time. We reinvestigate the 
case of Einstein gravity and point out additional features. Then we extend our results to include Gauss-Bonnet 
corrections. 

Our paper is organized as follows: we present the model in Section 2 and discuss the planar limit of our 
solutions in Section 3. We then present soliton solutions with scalar hair in Section 4 and "hairy" black holes 
in Section 5. We conclude and summarize in Section 6. 



2 The model 

In this paper, we are studying the formation of scalar hair on electrically charged black holes and solitons in 
(4 + l)-dimensional Anti-de Sitter space-time. The action reads : 

S = ^^ f d 5 x^=g~ (R - 2A + | (R mnkl Rmnkl - AR MN R MN + R 2 ) + 16^G£ mattor ) , (1) 

where A = —6/L 2 is the cosmological constant, a the Gauss-Bonnet coupling and M,N,K,L = 0,1,2,3,4. 
^matter denotes the matter Lagrangian : 

C matteT = --F MN F MN -(D M i')*D M il;-m 2 ^*^ , M,N = 0, 1,2,3,4 , (2) 

where Fmn = 9m An — On Am is the field strength tensor and DmtP = dMip — ieAMip is the covariant derivative, 
e and m 2 denote the electric charge and mass of the scalar field ip, respectively. 

The coupled gravity and matter field equations are obtained from the variation of the action with respect 
to the matter and metric fields, respectively, and read 

Gmn + Ag MN + -^Hmn = S7rGT MN , M,N = 0, 1, 2,3,4, (3) 



where Hmn is given by 

Hmn = 2(RmabcRn —ZRmanbR — ZRmaRn + RRmn) — -^9mn (R — ^RabR + RabcdR ) 

(4) 



and Tmn is the energy- momentum tensor 

, C^tnatter 

~dg~ 
We choose the following Ansatz for the metric : 



m r r> "•'-matter /-■> 

J-MN = 5M AT •'-matter — ^ a~MN ' ^ ' 



ds 2 = -f(r)a 2 (r)dt 2 + j^-dr 2 + ^dO| (6) 

where / and a are functions of r only. For the electromagnetic field and the scalar field we have [1] : 

A M dx M = </>(r)dt , 1> = il)(r) (7) 

such that the solutions possess only electric charge. 
The equations of motion then read : 

1 - / + 2r 2 /L 2 rf_ f 2e 2 ^ 2 + f(2m 2 a 2 ^ 2 + ft 2 ) + 2f 2 a 2 jj' 2 \ 

! ~ T r 2 + 2a(l-f) 1 2fa 2 \ r 2 + 2a(l - /)) ) (j 

, r 3 (e 2 ^ 2 + a 2 /V 2 ) 

a - 7 af 2 (r 2 + 2a(l-/)) W 

3 a'\ , e 2 V> 2 , N 

raj f 



3,/' a'N., /eV 

^ + 7 + aJ^"l7^ 2 / 



V/' - -(~ + ^-W'-(^-^W (ID 



where 7 = 167rG. Here and in the following the prime denotes the derivative with respect to r. These equations 
depend on the following independent constants: Newton's constant G, the cosmological constant A (or Anti-de 



Sitter radius L) and the charge e and mass m of the scalar field as well as on the Gauss-Bonnet coupling a. 
The system possesses two scaling symmetries: 



as well as 



Xr , t—*\t , L -¥ XL , e^f e/X , a -¥ X a (12) 



cf> ->• X(j> , V ->■ M> > e -> e/A i 7 -> 7 A 2 (13) 



which we can use to set L = 1 and 7 to some fixed value without loosing generality. 

Asymptotically, we want the space-time to be that of global AdS, i.e. we can choose a(r — >• 00) — » 1. Other 
choices of the asymptotic value of a(r) would simply correspond to a rescaling of the time coordinate. The 
matter fields on the other hand obey : 



.2 , . t\ $- , Vm- 



(r > 1) = ix - Q/r z , V(r > 1) = -T- + -T- (14) 



with 



2q 



A_=2-^4 + m 2 L 2 ff , A + =2 + ^4 + m 2 L 2 ff , £ c 2 ff = - _ -= == ~ £ 2 (l - Q /L 2 + Q(a 2 )) . (15) 

Note that the value of the Gauss-Bonnet coupling a is bounded from above : a < L 2 /4 where a = L 2 /4 is the 
Chern-Simons limit. In the following we will be interested in two cases: that of a tachyonic scalar field with 
to 2 = — 3/L 2 and that of a massless scalar field with m 2 = 0. Note that while for the tachyonic case, the fall-off 
of the scalar field will depend on a, it is independent of a in the massless case. 

The parameters /i, Q are the chemical potential and the electric charge, respectively. In the following we 
will choose ip- — 0. -0+ will correspond to the expectation value (O) of the operator O which in the context 
of the gauge theory-gravity duality is dual to the scalar field on the conformal boundary of AdS. Note that in 
comparison to the case of holographic superconductors (see e.g. [4] and reference therein) the space-time in this 
study possesses spherically symmetric sections for constant r and t. In other words, for r —> 00 our space-time 
corresponds to global AdS. 

The energy of the solution is given is terms of the coefficients that appear in the asymptotic fall-off of the 
metric functions which reads 

2 r 

f(r > 1) = 1 + -£- + ^ + • • • . a(T » 1) = 1 + °-\ + ■ ■ ■ ■ (16) 

h eff r r 

f% and C4 are constants that have to be determined numerically and that depend on the couplings in the model. 
The energy E then reads 



V 3 o,.r ,. h 



E= 8^G 3M ' M = -yVl-4a/£ 2 • (17) 

No explicit solutions to the system of equations (j5])- (fTTj) are known for ip(r) ^ 0. The solutions have to be 
constructed numerically. We have done this by employing a collocation method for boundary-value ordinary 
differential equations, equipped with an adaptive mesh selection procedure [37) . In the following, we will use 
the rescalings (|12p . (fT3"f and choose L = 1 and 7 to some fixed value. Note in particular that the scaling (IT31 
leads to e 2 — > e 2 /j. Due to convenience, we will choose 7 = ^ in the following. Hence, our e 2 will be related 
to the e 2 of [32] - which we denote e 2 in the following - by e 2 = ^e 2 . 

3 The planar limit 

In [36] it was pointed out for the a = limit and in 4-dimensional global AdS that the solutions can be 
connected to the planar counterparts in the limit of large charge. This is possible since the electromagnetic 
repulsion balances the gravitational attraction present in AdS space-time. As such the solutions can become 
comparable in size to the AdS radius L. On the AdS boundary these solutions look similar to the corresponding 



planar counterparts. We find that this is also true in 5-dimensional global AdS and including Gauss-Bonnet 
corrections and we will present our numerical results for soliton solutions in the follow section. Here, let us 
demonstrate the idea using an analytic and electrically charged black hole solution to the equations of motion 
for ip(r) = 0. This reads QH H3 H3 CSJ H2 

, a(r) = 1 , 0(r) = % - § , (18) 

where M and Q are arbitrary integrations constants that can be interpreted as the mass and the charge of the 
solution, respectively. In the limit a — > 0, the metric function f(r) becomes /(r) = 1 + fj — 73- + 4- and 
the corresponding solutions are Reissner-Nordstrom-Anti-de Sitter (RNAdS) black holes. First note that for 
M = Q = this solution corresponds to global 5-dimensional AdS with AdS radius L e g. Now let us rewrite the 
metric function f(r) as follows 





f(r)=r< -, + - l-\l- T , + — i -SP- (19) 



and adapt a rescaling to the conformal AdS patch that reads 

r -> Xr , i -> \~H , M -► A 4 M , Q -> A 3 Q , (20) 

where A is some arbitrary scaling parameter. For A — > 00 this tends to the planar solution with X 2 dfl 2 — > 
dx 2 + dy 2 + dz 2 . This corresponds to M and Q going to infinity. In addition we have the following rescalings 

(i-tXfi , i/>± -> A A± V± • (21) 

In other words, whenever we have a solution in global AdS (black hole or soliton) for which any of the parameters 
above tends to infinity we approach the planar limit. 

4 Solitons 

To find soliton solutions of the equations of motion, we have to fix appropriate boundary conditions at the 
origin r = 0. These read 

/(0) = 1 , 4/(0) = , V'(0) = , o'(0) = (22) 

and ensure that the metric and matter functions are regular at the origin. The solitons can then be characterized 
by the values of the matter and metric functions at the origin (f)(0), ip(0) 7 a(0) which depend on the choice of 
e 2 , Q, a and m 2 . 

For a = the existence of solitons in global AdS was discovered in j33j , where a perturbative approach was 
taken. In |32j it was shown that solitons can have arbitrarily large charge Q for large enough gauge coupling 
e 2 , while for small gauge coupling the solutions exhibit a spiraling behaviour towards a critical solution with 
finite charge and mass. In both cases the scalar field was chosen to be massless, i.e. m 2 = 0. The results of [32] 
were extended to the case m 2 = —3 in 35j. The qualitative results remain similar. In the following, we will 
reinvestigate both m 2 = and m 2 = — 3 and point out further details. 

4.1 m 2 = -3 

4.1.1 Charged solitons for a = 

The numerical results demonstrate that for fixed Q the solitons exist only for large enough values of e > e c , 
i.e. exist on an interval e G [e c (Q) : 00]. This seems natural when considering that small e 2 corresponds to 
large 7 and vice versa. Hence, solitons exist only if the gravitational coupling is small enough in comparison 
to the gauge coupling e 2 . An estimation of the critical value e c (Q) for Q £ [0 : 10] is presented by the solid 
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Figure 1: The value of the charge Q as function of the value of the gauge coupling e 2 at which a(0) = for 
soliton solutions with m 2 = and a = (dashed), m 2 = —3 and a — (solid black), m 2 = —3 and a = 0.1 
(solid red), respectively. The dotted-dashed red curve is just an upper estimate of the critical value e 2 since the 
numerical integration becomes very difficult here. Note that soliton solutions exist only to the right, below and 
above these curves. 



black line in Fig. [T] At the approach of e c (Q), our numerical results show that the value of the metric function 



a(r) at the origin, a(0) tends to zero, accordingly the metric becomes singular. In Fig 2(a) the quantity a(0) is 
plotted as function of e 2 for different values of Q. We find in particular e 2 « 1.04, 1.32, 1.36 for Q — 0.5, 1.0, 
2.0, respectively. This can also be seen in Fig. Q] The mass and the value "0(0) corresponding to these families 
of solutions are presented in Fig. 2(b) Both M and 0(0) remain finite for e — > e c . However we observe that 



V>"(0) becomes quite large and likely tends to infinity in this limit. 

These figures reveal that the dependence of a(0) on the coupling constant e at the approach of the critical 
value is not simple because there exist several branches of solutions, i.e. for a fixed value of e 2 more than one 
soliton exists and these solutions are distinguished by their value of a(0). Within our numerical investigations, 
we were able to produce up to three solutions with the same e and different values of a(0) but it might be that 
more oscillations occur while decreasing a(0) (we investigated solutions down to a(0) ~ 10 -5 ). Note that this 
behaviour is also present in the case m 2 = [32] . 

Profiles of solitons close to the critical limit are shown in FigJ3J More precisely, the metric and matter 
functions of the solutions with a(0) = 0.01 and a(0) = 0.001, respectively are shown for Q — 0.5 (corresponding 
to e 2 sa 1.043). The functions do not change much with a(0) apart on a small interval close to the origin. The 
metric function f(r) first decreases away from its value at the origin /(0) = 1. Then it reaches a minimum 
while staying nearly constant on a plateau, finally it increases to reach its asymptotic behaviour f(r) ~ r 2 . 

Our numerical results confirm that the solitons can be constructed for arbitrarily large values of e. In the 
limit e — > oo the metric functions tend uniformly to f(r) = 1 + r 2 and a(r) = 1. Note that the limit e 2 — > oo 
corresponds to a decoupling of the matter from the metric fields, i.e. 7 = (the so-called "probe limit" in 
the context of holographic superconductors). Hence, we have a fixed global AdS background and two coupled 
matter field equations for and ip (after suitable rescalings). 

For the purpose of later comparison of the pattern of the solutions for a = and q /0 (see below) we 
present the dependence of e 2 , of 0(0) and of ip+, respectively on a(0) in FigH] Here we have chosen Q = 0.5. 
We observe that all these quantities stay finite at the approach a(0) — > if a — (red curves). We will see in 
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Figure 2: The value of the metric function a(r) at the origin, a(0) as function of e 2 (left) and the value of the 
scalar held function ip(r) at the origin, ipi®) an d of the mass parameter M of the solution as function of e 2 
(right). Here Q = 0.5, 1.0, 2.0, respectively, a = and m 2 = —3. 



the following that this is different for a ^ 0. 



4.1.2 Charged Gauss-Bonnet solitons 

The numerical results strongly suggest that the soliton solutions which exist for a = get smoothly deformed 
when we gradually increase the Gauss-Bonnet coupling constant a. These branches of Gauss-Bonnet solutions 
exist up to the Chern-Simons limit a — L 2 /4 = 1/4. Figfl] suggests that soliton solutions exist for a larger 
parameter range in Q and e 2 when choosing a ^ 0. Note that solitons exist to the right, above and below the 
curves. Hence, the domain in the e 2 — Q— plane, in which soliton solutions do not exist is much smaller for 
a = 0.1 (red line) as compared to a — (black line). In the following, we will describe the properties of these 
charged Gauss-Bonnet solitons in more detail. 

In Fig 5(a) we give a(0) as a function of e 2 for different values of a and Q = 1. This is to understand how 
the critical behaviour in the limit a(0) — > changes when including Gauss-Bonnet corrections. We observe that 
increasing a from a = the value of e c first increases. Again, this results in the appearance of several branches 
of solutions (up to three for a = and a = 0.01). When a is large enough, only one branch of solutions exists, 
for our choice of parameters this happens at a > 0.05. Increasing a beyond a — 0.05 the critical value e c 
decreases strongly. Although this pattern was checked only for a few values of Q, we believe that it is generic. 

The masses corresponding to the various branches of Fig. 5(a) are given in Fig. |5(b)| It turns out that 
the mass of the soliton increases when increasing the Gauss-Bonnet parameter a and/or decreasing the gauge 
coupling e 2 such that the solution with the smallest e 2 and largest a has the biggest mass. Note also that for 
a = and a = 0.01, respectively the three branches of solutions mentioned above are shown in this figure. 
The masses of these branches do not differ strongly, this is why they cannot really be distinguished in the 
plot. However, we observe that for a fixed values of e 2 the solution on the first branch, i.e. the solution with 
the largest value of a(0) has the lowest mass, while the solution on the third branch, i.e. the solution with 
the smallest a(0) has the highest mass. This suggests that the solutions on the second and third branch are 
unstable, while the solutions on the first branch are stable. 

is qualitatively different from the 

0.5. 



case a 



Fig. 5(a) suggests further that the approach of the critical limit a(0) 

0. This can clearly be seen in FigQ] where e 2 , 0(0) and -0+ are given as function of a(0) for Q 
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Figure 3: The metric and matter functions of a soliton solution with Q = 0.5, e 2 rj e 2 , a = 0, m 2 
a(0) = 0.01 (red) and a(0) = 0.001 (black), respectively. 
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Figure 4: We show e 2 , ip+ and i/j(0), respectively as function of a(0) for Q = 0.5 for a = (red) and for a = 0.1 
(black). 




(a) a(0) 



(b) M 



Figure 5: The value of the metric function a(r) at the origin, a(0) as function of e 2 (left) and the value of the 
mass parameter M of the solution as function of e 2 (right) for several values of a. Here Q — 1.0 and m 2 = —3. 



As mentioned above, e 2 , 4>(0) and ip+ stay finite for a(0) — >• if a = 0. For a = 0.1, on the other hand, the 
decrease of a(0) results in a strong decrease of all three quantities. 

To understand this behaviour better, we have looked at the profiles of the metric and matter functions in this 
critical limit. In Fig. [5]we show the profiles of Gauss-Bonnet solitons with a — 0.1 and Q = 0.5 approaching the 
critical limit a(0) ->• 0. The two solutions correspond to a(0) = 0.01 (e 2 = 0.29) and a(0) = 0.001 (e 2 = 0.45), 
respectively. As can be clearly see in this figure, the phenomenon observed for a = (see Fig. [3]) exists also 
here and is in fact enhanced: we find that the metric function f(r) possesses a very pronounced local minimum 
which approaches zero for decreasing a(0). If the value of /(r) at the minimum reached zero before a(0) then the 
corresponding limiting solution would correspond to an extremal charged Gauss-Bonnet black hole. However, 
our numerical accuracy does not allow us to confirm this statement. At the moment that the value of a(0) 
reaches a value that is smaller than our numerical accuracy f(r) at the minimum has not yet reached zero. 
We find e.g. for this choice of parameters that for a(0) = 0.01 the minimum is located at r = r m « 0.317 
with a{r m ) w 0.54, f(r m ) = 0.06 and ip(r m ) = 1-7. For o(0) = 0.001 we find r m « 0.309 and a{r m ) « 0.49, 
f( r m) = 0.01, ip(r m ) — 1.46. So we cannot make a clear statement about the nature of the limiting solution 
yet. This is a difficult numerical question that we plan to investigate further in the future. However, we are 
sure that the limiting solution will not be a extremal Gauss-Bonnet black hole with scalar hair since it can be 
proven (see Appendix) that extremal Gauss-Bonnet black holes do not support scalar hair. The appearance of 
the minimum of f(r) is connected to a strong variation of e 2 when decreasing a(0). We further observe that 
a(r) remains roughly constant in the central region of the soliton and varies strongly in the region where f(r) 
attains its local minimum. The pronounced dependence of the electric potential on the value a(0) also contrasts 
with the a = case. 

The observation that for a fixed value of e 2 and of a several branches of solutions can exist leads to a new 
type of pattern when plotting a(0) as function of a for fixed values of e 2 and Q. This is shown in Fig. [7] for 
Q = 1 and three values of e 2 . For small gauge coupling (for our choice of parameters e 2 < 1.28) solitons do not 
exist for a = 0. However, if e 2 is large enough (here e 2 > 1.28) solitons exist for the full range of a € [0 : 0.25]. 

For e 2 = 1 the value a(0) is a monotonically increasing function of a with a(0) sw 0.1 at a — 0.25 and 
a(0) = at a critical value a c (e 2 ) « 0.03. For e 2 = 1.5, on the other hand, the value a(0) is a monotonically 
decreasing function of a with a(0) sw 0.49 for a — and a(0) s=s 0.16 for a = 0.25. For e 2 = 1.4 we observe 
that there is a main branch on which a(0) is monotonically decreasing with a with a(0) « 0.425 for a = 




Figure 6: The metric and matter functions of a soliton solution with Q = 0.5, a = 0.1, m 2 
(o(0) = 0.01) (red) and e 2 = 0.29 (o(0) = 0.001) (black), respectively. 
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Figure 7: The value of the metric function a(r) at the origin, a(0) as function of a for Q = 1, m 2 = —3 and 
several values of e 2 . 
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and a(0) w 0.15 for a = 0.25. However, there also exists a second branch that is disconnected from the first 
and has the form of a small loop that is present for small values of a and a(0). This is in fact a result of the 
complicated structure of the branches of solutions that exist for this choice of parameters. This small extra 
branch disappears when increasing e 2 . 



4.2 



m 







This case has been studied in detail in 32 for a = 0. It was found that for e < e 2 = 4p soliton solutions 
exist only for small values of Q and that at the maximal possible value of the charge Q the value of the metric 
function a(r) at the origin, a(0) tends to zero. For e 2 > 32/2 on the other hand solutions can exist for arbitrarily 
large Q with Q diverging for a(0) — > 0. In the following we will demonstrate that for e 2 < ^ there exists in 
fact a second branch of solutions at large Q which is separated from the first branch by a "forbidden band" 
of Q-values. The existence of a second branch of solutions was already noticed in 36 for solitons in global 
4-dimensional AdS. Note that our e 2 is related to e 2 by e 2 



-^e 2 such that 
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4.2.1 Charged solitons for a = 

The domain of existence of solutions in the case m 2 = is given in Fig[T] (black dashed line). As already 
mentioned previously solitons exist to the right, above and below this line which itself gives an estimation on 
the value of the charge Q in dependence on e 2 at which the value of a(r) at the origin vanishes, i.e. where 
a(0) = 0. For m 2 = and a = we find that hairy soliton solutions exist for Q £ [0 : oof if e 2 > e 2 - confirming 
the results stated in [32]. In agreement with J32] we find that e 2 « 2.4. 

For e 2 > e 2 > e 2 '* > on the other hand soliton solutions exist for Q £ [0 : Qf,i] and for [Q/.2 : oof. The 
difference between Q/,i and Q/,2 becomes larger when decreasing e 2 . At e 2 '* we find that Q/,2 — >• oo such that 



soliton solutions exist only for Q £ [0 : Q/,i]. 
a = 0. 



Our numerical results show that e 2 '* w 0.94 for m 2 — and 
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(b) M as function of i/>(0) 



Figure 8: The value of the charge Q as function of the value of the metric function a(r) at the origin, a(0) (a) 
as well as the mass M as function of the value of the scalar field function ip(r) at the origin, i/j(Q) (b) for soliton 
solutions with m 2 = 0, a = 0, e 2 = 2.0 < e 2 , e 2 = 2.37 < e 2 , e 2 = 2.42 > e 2 and e 2 = 2.5 > e 2 , respectively. 



In Fig 8(a) we show the value of Q as function of a(0) for e 2 = 2.0 < e 2 , e 2 = 2.37 < e 2 , e 2 = 2.42 > e 2 



and e = 2.5 > e 2 , respectively. As stated above, we find that for e < e 2 a second branch of solutions exists. 
For small e 2 this is separated from the small Q-branch by a "forbidden band" Q £ [Q/.i : Q/,2] in which 
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Figure 9: The value of the metric function a(r) at the origin, a(0) as function of e 2 for a massless scalar field 
and different values of Q and a. 

soliton solutions do not exist. Increasing e 2 these two branches approach each other such that at some stage 
they intersect for some given value of a(0). This is seen clearly for e 2 = 2.37. Increasing e 2 even further such 
that e 2 — > e 2 we find that the two branches start to merge and become one continuous curve for e 2 > e 2 . For 
e 2 > e 2 there are no solutions for a(0) = 0. In Fig |8(b)| we show the corresponding mass M of the solutions in 
dependence on tp(0). These curves are qualitatively similar to the curves in 4-dimensional AdS given in [36 . 
The M = Q = limit corresponds to global AdS with radius L e g. This corresponds to a(0) = 1 in the Q-plot 
and ip(0) = in the M-plot. Moreover, we observe that solutions with Q and M tending to infinity are possible. 
These are the planar solutions mentioned previously. Moreover, for ijj(0) becoming large (and equivalently a(0) 



tending to zero for 



< e; 



and a(0) tending to its critical, non- vanishing value for e 2 



:), the mass of the 



solutions tends to a constant. 

The critical behaviour at the approach a(0) — > was discussed in [35] and is qualitatively similar to the case 
to 2 = —3 discussed in this paper. In Fig|H]we show the dependence of a(0) on e 2 for different values of Q (black 
lines). 

4.2.2 Charged Gauss-Bonnet solitons 

We observe that the dependence of a(0) on e 2 looks qualitatively similar to the case m 2 = —3. This can be seen 
in Fig 1^1 where we show a(0) as function of e 2 for several values of a and of Q. Increasing a from a = we 
observe that the oscillations around the critical value e c are first amplified (see red curve in Fig|§]for a = 0.01 
and Q = 0.5). This results in the existence of several solutions with different values of a(0) for the same value 
of e 2 . E.g. we were able to construct up to three solutions for a = and a = 0.01, respectively. Increasing a 
further, only one branch of solutions exists. For our choice of parameters this happens at a > 0.05. Increasing 
a beyond a — 0.05 the critical value e c decreases strongly while increasing a. 

5 Black holes 

In order to find an explicit black hole solution of the equations of motion, we have to fix appropriate boundary 
conditions. In the following, we are interested in the formation of scalar hair on electrically charged black holes 
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CO 




Figure 10: The value of a(rh) as function of r^ for m 2 = — 3/L 2 = —3, e 2 = 1, Q = 1 and different values of 
the Gauss-Bonnet coupling a. 



with a regular horizon at r = r^ such that 

f{Th) = 

with a{rh) finite. In order for the matter fields to be regular at the horizon we need to impose: 



<f>(r h ) = , i>'(r h ) 



m 2 ^ (r 2 + 2afc) 



2r + 4r/L 2 - 7 r 3 (m 2 ^ 2 + (j>' 2 /(2a 2 )) 



(23) 



(24) 



Black hole solutions can therefore be characterized by the parameters r/,, Q, e 2 as well as the Gauss-Bonnet 
coupling a. 

The Hawking temperature of these black holes reads 



Tu 



while the entropy is 



1 

-in 



V-9 u 9 KL 9K9ttd L gtt 



1 

4ji 



f(r h )a(r h ) , #,£=1,2,3,4, 



S 



2G 



dark) 



(25) 



(26) 



The stability of black hole solutions for a = and m 2 = was studied in 32 . It was found that RNAdS 
black holes are never unstable to scalar condensation for e 2 < 3 and unstable for all charges Q for e 2 > ^. In 
the intermediate interval for 3 < e 2 < ^ RNAdS black holes become unstable for sufficiently large charge Q, 
i.e. for Q > Q c (e 2 ). Moreover, it was observed that for large gauge coupling and small charges the solutions 
exist all the way down to vanishing horizon. The limiting solutions are the soliton solutions mentioned above. 
On the other hand for large charge the limiting solution is a singular solution with vanishing temperature and 
finite entropy, which is not a regular extremal black hole |34) . Similar results were obtained for a tachyonic 
scalar field with m 2 — —3 in 35 . Here we reinvestigate both cases and point out further details. 



5.1 m z = -3 

5.1.1 Hairy charged black holes for a 







We have first studied the interval of ru on which black holes exist for Q = 1 = e 2 . This is shown in FigflOl 
where we give a(r/j) as function of r^. At the maximal value of r/i, rh,max the solution tends to a black hole 
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Figure 11: The matter and metric functions of a black hole solution with r^ ~ 0.01 and the corresponding 
soliton solution for Q — 1, e 2 = 2. 

solution without scalar hair, i.e. the RNAdS solutions with a(r) = 1 and ip(r) = 0. The value rh,max depends 
on e and Q, e.g. for Q = 1 we have rh, ma ,x = 0.56, 0.59, 0.62, respectively for e 2 = 1.0, 1.5, 2.0. 

At the minimal value of rh, r^^mm the solutions either tend to a singular black hole solution with a(rh, m in) = 
and correspondingly temperature Tr = or to the soliton solutions mentioned early in this paper. For fixed 
Q this depends on the choice of e 2 . This can been see in FigfTUl and FigfTTJ The curve of a(0) tends to zero for 
r/j > when choosing Q = l = e 2 ,a = (see Fig flO]) . Increasing e 2 this changes. For e 2 large enough the black 
holes exist all the way down to r^ = 0, where they merge with the soliton solutions discussed above. Generally, 
we can say that for a fixed Q the black holes tend to the soliton solutions for e 2 > e 2 (i.e. for values of e 2 for 
which soliton solutions exist), while they tend to the singular black hole solutions for e 2 < e 2 . For Q = 1 we 
have found (see discussion on solitons above) that e 2 « 1.32. Hence we would expect that the black holes merge 
with the solitons for e 2 > 1.32. This is indeed what our numerics shows and is given for Q = 1 and e 2 = 2 in 
Fig. [TT] Here, we plot the matter and metric functions corresponding to a black hole with small ru = 0.01. 
Clearly, the black hole solution approaches the soliton for rh — > 0. While a(r) and ip(r) are nearly identical 
on the full r-interval, the behaviour of f(r) and <p(r) close to the origin differs between the black hole and the 
soliton case. This is not surprising, since the boundary conditions for black holes are f(rh) = 0, 4>{rh) = 0, 
while for solitons we choose /(0) = 1 and 4>(0) = c/)q > 0. Fig fTT] also suggests that f'(rh) and 4>'(rh) tend to 
infinity for r^ — >• 0. This is shown in Fig ll2l where we give the values of tfj(rh), f'(vh) and (f>'(rh), respectively 
as function of rh for e 2 = 2 and Q — 1 (black lines). Clearly, both /'(r^) and 4>'(rh) tend to infinity. The value 
tpifh) reaches a finite value for r^, — > 0, e.g. ip(rh — > 0) — > 1.4 for e 2 = 2, Q = 1. 

We have also studied the behaviour of the solutions for fixed r^ and varying e 2 . We find that the black 
holes exist for large enough values of e 2 . In the limit e 2 — > oo the scalar field function slowly approaches the 
null function. However, we observe a new phenomenon: we find that the metric function f{r) is no longer 
monotonically increasing and develops a local minimum at r m where r^ > r m > oo for sufficiently small e . 
This is shown in Fig. [13] for rh = 0.525, Q — 1 and e 2 = 0.6975. Decreasing e 2 the value of the local minimum 
of f(r) decreases further and the solution evolves into a configuration with different features for r G [rh, r rn ] and 
r G [r m ,oo], respectively. For r G [r^,r m ] the scalar field function ip(r) is clearly non- vanishing and the value 
a(rh) becomes very small - in fact of order 10~ 6 . This renders the numerical construction quite involved and it 
becomes unreliable at some stage because the value a(rh) gets smaller than the accuracy of our numerical code. 
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Figure 12: The values of ip( r h), f'( r h) and <f>'(rh), respectively as function of r/, for a = (black) and a = 0.01 
(red). Here e 2 = 2, Q = 1. 

On [r m , oo] we have a(r) ~ 1 and ip(r) << 1. 

The occurrence of this phenomenon depends crucially on the fact that a RNAdS black hole exists for the 
chosen values of Q and r/,. RNAdS black holes with horizon at r — rn indeed exist for small Q and become 
extremal for Q = Q m with Q 2 m = r^(l + rj^/j. When considering hairy black holes with Q and rt we find the 
following when decreasing e: (a) if a RNAdS solution with the same Q and r/i exists then the hairy black holes 
tend to the corresponding RNAdS black hole for a finite value of e, (b) if a RNAdS black with the same Q and 
Th does not exist an additional local minimum develops and the above described phenomenon appears. 



5.1.2 Hairy charged Gauss-Bonnet black holes 

Our numerical results suggest that the black holes which exist for a = get continuously deformed for < a < 
L 2 /4 = 1/4. In the following, we want to discuss in which way these charged Gauss-Bonnet black holes differ 
from the charged black hole solutions in the a = limit. 

We have first studied the interval of the horizon radius r/j, [rh,min '■ fh,max] on which the black hole solutions 
exist for fixed Q, e 2 and a. At the maximal value of r/i, rh, m ax the solution tends to a black hole solution 



without scalar hair, in which case ip(r) = and the remaining functions are given by (|18[) . 

At the minimal value of r^, rh, m in the solutions either tend to a singular black hole solution with a(rh, m in) = 
or to a singular solution with /(r^) — > for r^ — > 0, while a(rh) stays finite. We will demonstrate in the 
following that the limiting solution for ru — > is not the Gauss-Bonnet soliton solution discussed above. 

In FigfTU]we give a(rh) as function of rn for different values of a for Q = 1 and e 2 = 1. We observe that 
Th,max increases slightly with a. Moreover, for Q = e 2 = 1 there is no soliton solution for a = and the black 
hole solutions end at a singular black hole solution. This however changes when increasing a. For a sufficiently 
large (for our choice of parameters for a > 0.01) there exist black hole solutions for r/, €]0 : rh,max\- This is 
clearly seen in FigflOl 

One important observation is that the limit r/,. — > is different as compared to the a = case. This is 
shown in FigfMl where we give the matter and metric functions of a Gauss-Bonnet black hole with r^ = 0.01 
and of the corresponding soliton solution for e 2 = Q = 1 and a — 0.24. Rather then tending to the soliton, 
the black hole approaches a singular configuration for 77,. — > 0. To strengthen this claim we show f'(rh), 4>'( r h) 
and ip(rh), respectively, as function of r% in Fig. Q21 For a/0we find that the values f'(rh), </>'(f/i) smoothly 
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Figure 13: The metric and matter functions of a hairy RNAdS solution with 77, = 0.525, Q — 1 and e 2 = 0.6975 
are given. 



reach zero for r% —> while ip(rh) diverges. As a consequence, the Gauss-Bonnet black holes do not approach 
the Gauss-Bonnet solitons discussed above in the limit ru — > 0, which have /(0) = 1 and 0(0) = <fio > 0. This 
is completely different then in the a = case where f'{rh) and 4>'(rh) both tend to infinity for r^ — > 0. 

The difference is also seen when comparing the entropy, temperature and mass of the limiting solutions. For 
a = there are two different limits: if the black hole approaches the soliton, the limiting solution has infinite 
temperature, while the entropy tends to zero. On the other hand, if the black hole tends to a singular black 
hole solution the temperature tends to zero, while the entropy stays finite. For a > the temperature also 
seems to tend to zero in the limiting case. Independent of the choice of a the mass remains finite. 

We also observe that ip + remains finite and positive in the limit r^ — > 0. The effect of a on the quantities 



given in Fig. [12] is shown in Fig 15(a) Complementary to the figures above, the dependence on a of the same 



parameters characterizing the solutions is shown in Fig. 15(b) 



The behaviour of the black hole occurring while the coupling constant e 2 decreases seems to persist for 
a > 0. This is illustrated in Fig. 1161 The numerical results suggest that a local minimum of the metric function 
f(r) forms that seems to tend to zero. Our numerical investigation doesn't allow as to make a clear statement 
about the limiting solution. However, given the prove in the Appendix, we are sure that this is not an extremal 
Gauss-Bonnet black hole with scalar hair. 



5.2 

5.2.1 



m 







Hairy charged black holes for a = 



We have also studied the dependence of our results on the mass of the scalar field. Our results for m 2 = 
(together with those for m 2 = —3) are shown in FigfTT] It is clear from this figure that for a — the results do 
not depend on m 2 . 

5.2.2 Hairy charged Gauss-Bonnet black holes 

In Fig ll7l we also show the dependence of a(rh) on r^ for a = 0.24. Comparing this curve with the corresponding 
curve for m 2 = —3, we find that rh tmax increases with decreasing m . At some intermediate r^ the curves for 
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Figure 14: The matter and metric functions of the Gauss-Bonnet black hole with rt = 0.01 as well as those of 
the corresponding soliton with e 2 = Q = 1 and a = 0.24. 

m 2 = and m 2 = —3 intersect such that for smaller and fixed r/j the value of a(rh) increases with decreasing 
m 2 . In contrast to the a = case there is hence a dependence on m 2 . 

6 Summary and Conclusions 

We have studied the instability of static and charged solitons and black holes with respect to the condensation 
of a massless and tachyonic scalar field, respectively. We have re-investigated the case of Einstein gravity and 
have extended the results to include Gauss-Bonnet gravity terms. In both cases we have addressed the problem 
of existence of hairy solutions by solving the full system of equations numerically, i.e. our solutions are therefore 
not perturbative solutions. 

For the Einstein case, we find new features as compared to those reported in the perturbative limit in [32] 
and propose the following scenario for the domain of existence of hairy black holes. For fixed Q there exists an 
rh ex which corresponds to the horizon radius of the extremal RNAdS solution. Hence, RNAdS solutions exist 
only for r^ > r^ ex . For r^ — soliton solutions exist only for e 2 > e 2 . We hence find 

• for fixed (and small) Q and e 2 > e 2 the black holes tend to the corresponding soliton solution with the 
same Q when decreasing r^ 

• for fixed Q and e 2 < e 2 the black holes tend to a singular solution when decreasing rt 



for fixed Q and r^ > r/,,, 
decreasing e 2 

for fixed Q and r^ < rh,e 
metric function /(r) . 



the black holes tend to the RNAdS solutions with the same charge Q when 



the black holes tend to a configuration that possesses a local minimum of the 



For a/0 the scenario is similar except for the fact that the Gauss-Bonnet black holes never tend to the 
corresponding soliton solution when decreasing r^ - even when the soliton exists for the given values of Q and 
e 2 . It would be interesting to check whether these new features in the Gauss-Bonnet case could be recovered 
within a perturbative approach comparable to that used in the Einstein case in [32] . 



17 









i \ black: « = o.oi 
i \ — — . red : a -0.10 

l\ A ^ 






♦ <r h W4 












K / \ 










|N \ 


. 


I /\ \ 




1 / \ \ 




\i x \ --. 


- 


( ^<>^ 


/ v ^ ^^- ~~>^ 












i / ~^xr~~~^ \ 




/ / __— — ^~"^ "~^^> 








/>^ \ 









100 - 




~"~"~-~\^^ red :r h = 0.1O 


10 - 

1 - 

0.1 - 




~~~-- 


^- ~~~~' i \^\ 


0.01 - 




t'frj 

f(r h ) 

f'frh) 




0.001 - 






| \ 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



(a) a fixed 



(b) rj, fixed 



Figure 15: The values of ip(rh), f'{fh) and 4>'{ r h) as function of rn for two different values of a (left) and as 
function of a for different values of ru (right). Here e 2 — 2, Q = 1. 

Furthermore, none of the limiting black hole solutions can be a regular extremal black hole with scalar 
hair. Here we have extended the results of 011 and have shown that the near-horizon geometry of extremal 
Gauss-Bonnet black holes does not support massless or tachyonic scalar hair. 

In extension to [32] we observe that for intermediate values of the gauge coupling hairy solitons exist if the 
charge Q is small enough or large enough and that a "forbidden band" of charges is present in which hairy 
solitons do not exist. 
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7 Appendix: A No-hair theorem for extremal Gauss-Bonnet black 
holes 

In [32] it was argued that extremal Reissner-Nordstrom-AdS black holes cannot carry scalar hair. Here, we will 
investigate the case including Gauss-Bonnet corrections. 

In the following, we want to discuss the equations in the near-horizon case. It was shown in [35J that the 
near-horizon geometry is AdS2 X S 3 . We assume the metric to have the following form 

ds 2 = vi (-p 2 dr 2 + \dp 2 \ + v 2 (dip 2 + sin 2 ip (d9 2 + sin 2 6dip 2 )) , (27) 

where v\ and v 2 are positive constants. 
We find that the equations then read 

3 . /10' 2 p 2 , 2 e 2 (/>V 2 2 \ 
A = 8ttG ( -^s- + — 4' 2 + 2 — + m V . (28) 

V2 \2 Vf Vl p z Vi J 

± + A = 8 „ G (-_^ + £> + fW _,,, 

V2 \ I v{ Vl p A V\ 



+ A = 8ttG ---4- + ^-4' 2 + I r - mV » (29) 



1 . 1 2o ^A0' 2 P 2 ,,2 e 2 ^ 2 ^ 2 , ,,\ 
+ A + — + = 8ttG --V - ^' 2 + ; " - to 2 V> 2 • (30) 

U 2 Vl ViV 2 \2 Vi Vi p z Vi J 

Note that H TT and H pp are vanishing and hence the equations are the same as for a = (compare to [341 ) . 
Combination of (QH)) and ([29]) yields 

/p 2 /2 e 2 2 ^ 2 \ , s 

= 16ttG — ^' 2 + ^H- ■ (31) 

We can hence draw the same conclusion as in the a = case, namely that since both terms on the rhs of ([3TT) 
are positive they must vanish and hence tp' = and <p 2 ip 2 = in the near-horizon geometry. Since — — A > 
and to 2 < 0, the case = is ruled out from ([28]). hence the only possibility is tp = in the near- horizon 
geometry and we conclude that extremal Gauss-Bonnet black holes can not carry massless or tachyonic scalar 
hair. Furthermore ([30]) is always fulfilled in this case. To see this note that vi is equal to the square of the AdS2 
radius, i.e. V\ = R 2 , which was computed to be given by R 2 — (r 2 + 2a)/ (4 + 12r 2 /L 2 ) [35J [55]. Moreover, 
v 2 = r h> where 77, is the horizon radius. Inserting these expressions into ([50]) it is easy to see that the equation 
is fulfilled and doesn't put further constraints on the parameters. 
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